Cryogenic small-signal conversion with relaxation oscillations in Josephson junctions 
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Broadband detection of small electronic signals from cryogenic devices, with the option of simple 
implementation for multiplexing, is often a highly desired, although non-trivial task. We inves- 
tigate and demonstrate a small-signal analog-to-frequency conversion system based on relaxation 
oscillations in a single Josephson junction. Operation and stability conditions are derived, with 
special emphasis on noise analysis, showing the dominant noise sources to originate from fluctuation 
processes in the junction. At optimum conditions the circuit is found to deliver excellent noise 
performance over a broad dynamic range. Our simple models apply within the regime of classical 
Josephson junction and circuit dynamics, which we confirm by experimental results. A discussion 
on possible applications includes a measurement of the response to a cryogenic radiation detector. 
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I. INTRODUCTION 

Cryogenic devices are widely used in a broad range 
of applications like radiation detection, quantum crypto- 
graphy, charge manipulation on the single-electron level, 
quantum Hall effect or in basic studies of mesoscopic 
transport. Measurement of the electronic properties of 
such devices usually requires sophisticated readout elec- 
tronics. Detection schemes where the samples at cryo- 
genic temperatures are remotely connected to room tem- 
perature electronics generally face the problem of re- 
duced frequency bandwidth due to the impedance of long 
readout lines. In addition, the risk of noise pickup on 
the lines is intrinsically increased. Alternatively, sig- 
nal readout relatively close to the sample can be ac- 
complished with 'Superconducting Quantum Interference 
Device' (SQUID) amplifiers, which perform very success- 
fully in many cases but require a delicate setup (shield- 
ing) and are usually constrained to commercially avail- 
able systems. Amplification or impedance transforma- 
tion on-chip or very close to the sample is also possible 
with a 'High Electron-Mobility Transistor' (HEMT), the 
dissipation of which may, however, quickly reach an un- 
acceptable level. 

In recent years it has been realized that probing the 
electronic transport in a cryogenic device with a radio- 
frequency (RF) signali*2*MiSi& may have considerable ad- 
vantages compared to direct signal readout, mainly due 
to a substantial extension of the bandwidth. In those 
schemes the power of the reflected (or transmitted) RF 
signal from a properly tuned tank circuit is related to 
the electronic state of the device under test. The circuit 
needs to be carefully designed to minimize back-action on 
the cryogenic sample. Operation at microwave frequen- 
cies also naturally opens a potential way for frequency- 
domain multiplexing^^ 

A promising readout scheme, which we present in this 
paper, consists of an on-chip analog-to-frequency con- 
verter, delivering a frequency signal of large amplitude 
which is easily demodulated with standard room temper- 
ature (phase-locked loop) electronics. It has the advan- 
tages of both the direct signal readout close to the sample 
and a large frequency bandwidth. Particularly, it is much 



easier to accurately analyze a frequency signal than to 
transmit a low-level analog signal through long readout 
lines and amplify it with room temperature equipment 
which typically shows inferior performance in terms of 
noise with increasing temperature. Our low-noise con- 
verter circuit is based on a hysteretic Josephson junction 
exhibiting relaxation oscillations. Related ideas using 
relaxation oscillations in Josephson junction were pro- 
posed for thermometry? or direct radiation detection^ 
both relying on the temperature dependence of quasipar- 
ticle population in the gap singularity peak of asymmet- 
ric junctions, and for the (double) relaxation oscillation 
SQUID;ii*i2iiii which is investigated and used as a mag- 
netometer. In our case the circuit converts an analog 
current signal into a frequency with acceptable linearity 
over a broad operation bias and dynamic range. 

In Sec. [n] we review the basic principle of a relaxation 
oscillation circuit and derive conditions for stable opera- 
tion. Results from the model are illustrated with experi- 
mental data. A thorough noise analysis with implications 
for the circuit's readout resolution is given in Sec. II I II 
An optimized low-noise configuration with numerical es- 
timates is considered in Sec. IIVI followed by a discussion 
on possible applications in Sec. [V] As an example we 
demonstrate the readout of a cryogenic radiation detec- 
tor. The paper concludes with Sec. IVII 



II. PRINCIPLE OF OPERATION 

We assume a Josephson device with normal resis- 
tance R n , critical current I c , junction capacitance Cj 
and superconducting energy gaps Aj.,A2 , where < 
|Ai — A2I <C Ai + A2 = eV g . It shall be connected in 
series with an inductance L and both shunted with a re- 
sistor R s , as shown schematically in Fig.^ The circuit is 
eventually current biased by a large resistor R\, ^> R s . A 
Josephson junction with a non- vanishing difference of the 
energy gaps shows a region of negative differential resis- 
tance in the current-voltage (/V") characteristics. Volt- 
age biasing the junction in that region, where R s acts 
as voltage source with R s <C R n , its operation is po- 
tentially unstable and the circuit can undergo relaxation 
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FIG. 1: Left: Circuit diagram of the relaxation oscillator. 
Right: Schematic IV characteristics of a Josephson junction 
with an (exaggerated) region of negative differential resis- 
tance. The four partial processes of an oscillation cycle are 
labeled A, B, C and D. 

oscillations>i2ii£*i& A relaxation oscillation cycle, which is 
displayed in the IV diagram of Fig. can be separated 
into four phases: 

• (A) Initially, when V s = R s Ib > R S I C is turned 
on, the Josephson junction is essentially a short 
(supercurrent branch) and the current through L 
increases with a time constant t sc = L/R s towards 
a value J max = V S /R S like I sc (t) = J max (l -e~* /rsc ) 
until reaching I c within a time 

r A =-r sc ln(l-^). (1) 

• (B) Because the junction was current biased during 
phase (A) via a high-impedance L it switches now 
to the quasiparticle branch by developing a voltage 
across Cj until it is charged to V g . The inductance 
holds the current constant if Cj is small enough, 
which means that in order to observe the "full 
swing" of the voltage oscillations, the inductive en- 
ergy / c 2 L/2 and the energy from the bias voltage 
V^Cj/2 must be sufficient to provide the charge on 
Cj with V g , i.e. I 2 C L + V?Cj > (V g - V s ) 2 Cj must 
be fulfilled. For the case of interest where V s -C V g 
this requirement is particularly true if 

°i « W • ( 2 ) 

9 

Another requirement is undercritical damping of 
the R s LCj circuit with L/(CjR 2 s ) > 1. How- 
ever, because R s <C V g /I c , comparison with Eq. 
shows that the undercritical damping condition is 
already implied by @- The voltage switching time 
t b is on the order of CjV g /I c which is negligibly 
short for Cj satisfying J5J. 

• (C) Similarly to phase (A) the current on the quasi- 
particle branch decays with r qp = L/R qp (where 
R qp is the corresponding resistance in that region 



of the IV characteristics including the shunt R s 
in series) from I c to I m i n = (V s — V g )/R qp like 
I qp (t) = (I c - Imin)e~ t/Tqp + Imin until reaching zero 
(or a local minimum close to zero) within time 

• (D) The capacitor is discharged again to zero volt- 
age according to the conditions in phase (B), but 
with a subtle difference regarding final locking to 
the zero- voltage state, as discussed in the noise sec- 
tion cnsi 

Neglecting the short voltage switching times of phases 
(B) and (D), the relaxation oscillation period is given 
by (t a + t c ). However, when biasing a Josephson junc- 
tion at V s <C V g , the oscillation dynamics are dominated 
by the process in phase (A) with t a 3> r c 0{r A y-). 
Therefore, the relaxation oscillation frequency is essen- 
tially given by 

/r » T -_-* h -(l-I). (4) 

where a — Ib/I c is the reduced bias current. A series 
expansion for R S I C <C V s <C V g yields 

'--£H-°(=)} < 5 > 

« j^. (a»l) (6) 

These equations describe an almost linear analog-to- 
frequency converter. The same result follows from an 
expansion around R s — *■ 0, relevant for the readout of a 
variable resistance device in place of R s ■ The current-to- 
frequency conversion factor is 

Figure [21 shows experimental relaxation oscillation 
data. The amplitude of V out corresponds to the gap volt- 
age V g , and the dynamics follow the model predictions. 
In the operation range R S I C < V s < 0.1 ■ V g the effective 
oscillation frequency (t a +t b )~ 1 deviates from linear- 
ity by less than 10% (a larger operation range can also 
be chosen with an easy subsequent linearization of the 
results according to circuit calibration). Linear extrapo- 
lation of t^ 1 towards V s = yields a frequency offset in 
agreement with Eq. JSJ. 

Relaxation oscillations in Josephson junctions can be 
analyzed in terms of subharmonics of the Josephson fre- 
quency. This implies that the number of Josephson 
oscillations per relaxation oscillation cycle be much larger 
than unity in order to prevent significant harmonic phase 
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FIG. 3: Measured relaxation oscillation frequencies versus 
bias voltage for different I c due to application of magnetic 
field. The values to the right-hand side of the curves are the 
effective critical currents in (fJ,A) as obtained from fitting the 
theory to the experimental data. Device and circuit parame- 
ters are the same as in Fig. [5] 



FIG. 2: Relaxation oscillation time-traces, plots (a) and (b), 
measured at bias voltages differing by a factor of two. Plot (c) 
shows the inverse of measured time constants t a , t c and of 
their sum as a function of V s . The two dots in (c) correspond 
to the signal traces (a) and (b). Data were taken at 80 mK 
with an Aluminum based Josephson junction with nominally 
I c — 58.3 fiA and a circuit consisting of L = 280 nH, R s — 
91 mil 



locking of the two oscillating processes. That sets a con- 
straint on the frequency response and we can write 



fr = 



Vs 



(n<j, > 1) 



and 



I C L = 0o"0 , (a, n$ > 1), 



(8) 



(9) 



where <po is the magnetic flux quantum. This argument 
is in line with the requirement 



0L = ^IcL » 1 
00 



as stated in the literature^ 

Modulation of 7 C by application of a magnetic field 
parallel to the Josephson junction results in a variation of 
the relaxation oscillation frequency according to Eq. JSJl . 
This offers a convenient way to tune the circuit's dynamic 
properties as well as to extend the operation range to 
lower 7f,. Figure[31shows measurements of f r as a function 
of V s for different 7 C . In order to take an 7 C modulation 
into account we denote K = I® / I c as the factor by which 
the nominal value ij? may be suppressed. In the limit 
T — > and equal superconductors with gap A the zero- 
field critical current*^ is 7° = 7rA/(2ei?„). Corrections 



due to small gap differences are safely neglected for our 
analysis and we can write 



7 r 7L 



4k 



(11) 



The dimcnsionless capacitance parameter /3 C , which is 
a measure for the damping strength of the junction, is 
given by 



(12) 



where uj p = (2eI c /hCj) 1 ^ 2 is the Josephson plasma fre- 
quency. For a weakly damped and hysteretic Josephson 
junction we should choose f3 c larger than unity. The in- 
equality (J2J in terms of Eqs. (|11I12I1 can be rewritten as 



(13) 



(10) and substituting I C L from © yields 



> (3 C 



(14) 



This result is consistent with the condition 3> 1 as we 
required for Eqs. (|8I9|I . 

Finally, a comparison between (|ll)fl and (|13|) shows that 
the latter is the more stringent of both conditions by the 
factor I3 c (Ak/tt) 2 > 1. Consequently, the very minimum 
of I C L is determined by relation (| 1 31) . which constitutes, 
together with (|14|) and j3 c > 1, the relevant conditions for 
proper observation of relaxation oscillations and which 
should help to choose appropriate circuit components. 
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III. NOISE AND RESOLUTION 
A. General 

In this section we list the significant current noise 
sources referred to the circuit input (i.e. at Spe- 
cial attention is paid to experimental mean fluctuations 
of the relaxation oscillation periods, which are denoted 
by 5r r . Assuming an analog signal which requires a band- 
width fb w in order to resolve its dynamics in time (i.e. 
a sampling time period f^), we measure N = f r /fbw 
oscillations per sampled signal. The relative accuracy of 
a measurement improves with N as 



Sr r 1 




(15) 



Because bias current fluctuations 81b are linear to fre- 
quency fluctuations Sf r according to Eq. J7J), we can also 
conclude in first order that 5r r /T r = Sh/h- This yields 
an expression for the rms current noise of the signal sam- 
pled at f bw : 



(th) fb ,„ = (SIi 




(16) 



Consequently, in the case of random and uncorrelated Sr r 
fluctuations we observe a white current noise spectrum 
with a density 



3b 



{SI b 



h Sr r 

y/fr T r 



(17) 



apparently independent of fb w . Because phase (A) of a 
relaxation cycle dominates the timing, we expect fluctua- 
tions in the critical current SI C oc Sr r to be a major origin 
of St t noise. 



B. Flicker 1// noise in the critical current 

The critical current of Josephson junctions can fluctu- 
ate due to stochastic charge trapping at defect sites in 
(or close to) the barrier, which are known as "two-level 
fluctuators" . A sufficiently large ensemble of such fluc- 
tuators generates a 1// spectrum, with significant con- 
tribution only at low frequencies. According to empirical 
model o 1 ?' 2 ? the critical current noise density jtif due to 
1/f flicker noise can be described by 



Jtif 



ijT 2 
AS 



A ■ 



(18) 



where A is the junction area and A w 8.16 x 10~ 24 m 2 /K 2 
is an average value obtained from collecting data over 
a wide range of different junction parameters^ Scaling 
with T 2 was confirmed 2 ^ for temperatures down to 90 mK 
(although the authors 25 found a higher noise level in their 
devices with A « 39 x 10~ 24 m 2 /K 2 ). 



C. Critical current statistics from thermal 
activation 

Escape from the zero- voltage state of a Josephson junc- 
tion due to thermal activation is a well-known and widely 
studied phenomenon. It can be treated for a large vari- 
ety of junction types and external conditions. For our 
noise analysis we can restrict ourselves to the simple 
"transition-state" model 2 ^ 2 ^ where a particle inside a 
well is thermally excited above the relative barrier poten- 
tial and irreversibly leaves the bound state. This model is 
appropriate for underdamped Josephson junctions and is 
justified by typical device parameters and experimental 
conditions as given in the numerical examples (sections 
HHandEI). Particularly, we assume intermediate opera- 
tion temperatures satisfying 



7 



kT 
~E~j 



« 1, 



(19) 



where Ej = hl c /2e denotes a Josephson coupling energy, 
sufficient to suppress the probability of retrapping from 
the running statef 2 ^* 2 ^ 2 ^ but at the same time not too low 
to prevent macroscopic quantum tunneling effects. Ac- 
cording to the model there is a nonvanishing probability 
for transitions from the superconducting to the resistive 
state at current values I m < I c . The lifetime of the zero- 
voltage state in a Josephson junction as a function of the 
reduced current i — I/I c can be expressed by 26 - 27 



Zir 



(20) 



where ui a = u> p (l — i 2 ) 1 / 4 is the "attempt frequency" of 
the particle in the well and Uo = 2Ej(y/l — i 2 — i arccosi) 
is the relative potential height of the next barrier in the 
Josephson junction "washboard" potential. 38 The prob- 
ability P(t) for the junction to have switched from the 
superconducting to the resistive state before time t is 27 



P(t) = l-expl- J T^I{t') dt'j 



(21) 



By assuming small fluctuations compared to I c and using 
approximations in the limit e = 1 — i 1, we can solve 
the integral by a similar approach to Ref. |23 which yields 



P( e ) = l - exp 



w P 7 



47r/ r (2 e )V4 



exp 



2(2e 



,3/2 



3" 



(22) 

The mean value (I m ) of the observed critical current and 
its standard deviation (SI m ) are calculated in the Ap- 
pendix and are found to be 



and 



(SI„ 



1 - - (7 In 17) 



v 2/3 



2/3 



(ln^/a 



(23) 



(24) 



5 



where 



E. Noise conclusions 



V 



(zTT/r) (2) 



(25) 



Hence, (SI m ) is essentially a function of 7 2 / 3 , with a weak 
dependence on (f r /uj p ). The approximations used for 
derivation of I|23I24|) are appropriate for 77 ^ 1. Simi- 
lar results were obtained in Ref. 29. The current noise 
density j m at the circuit input is, according to Eq. lfT7|) : 



v 2 /3 



Jm 



(26) 



To verify our results and to compare with other 
models2S*2i of different formalism or treating different 
ranges of damping strength, we evaluated numerically 
the probabilities P(i), the transition current distribu- 
tions p(i) = dP(i)/di, and analyzed them with respect 
to shape, expectation value mi and width a m . We found 
that, within the range of allowed and reasonable model 
parameters, only the mean values m\ differed quanti- 
tatively for different models, as should be expected for 
different initial conditions and excitation forms. How- 
ever, there were insignificant differences in the distribu- 
tion shapes and particularly of their widths a m . There- 
fore, Eq.[^]can be considered a good estimate for critical 
current fluctuations due to thermally activated escape, 
applicable over a wide range of (3 C > 1. 



D. Other noise sources 

Thermal current noise from ohmic resistors is domi- 
nated by the shunt R s and corresponds to the standard 
Johnson noise j 2 = 4kT/R s . The voltage noise gener- 
ated by R s and seen by the junction is, due to f r cx V s , 
not amplified and therefore equivalent to j s at the circuit 
input. 

The real part of a good inductance L vanishes. There- 
fore, L can safely be considered as a "cold resistor" with- 
out thermal noise contribution. Pickup of external mag- 
netic noise can be shielded and becomes negligible for 
small coils. 

Because the Josephson current is a property of the 
ground state of the junction, it does not fluctuate. Hence, 
shot noise in Josephson junctions is only due to the quasi- 
particle current. The relaxation oscillations within our 
concept are dominantly determined by processes with the 
junction in the superconducting state. Therefore, shot 
noise by itself should be negligible in our case. 

However, as a consequence of the random nature of 
the junction phases in the quasiparticle tunneling regime, 
the locking to the zero-voltage stated at the end of an 
oscillation cycle occurs within a time spread^ on the 
order of 8t z rj 2%wLCj. This results in an input current 
noise density 



h St z 



(27) 



Combination of all noise sources derived above (and 
assumed to be uncorrelated) yields a total circuit input 
current noise density j b with j\ = j 2 l{ + j 2 ^ + j 2 + j 2 . 
We make substitutions with respect to a notation of % 
in terms of the primary circuit and operation parameters 
L, R s , R n , I c , a and T: 



. 2 {I c R n Tf , alTT^L 
Jb = c i ; ^c 2 - 



/ 



c 3 



T 
R~ s 



c 4 - 



a 3 I 2 R s 



Rn 



(28) 

where the constant coefficients Cj are orthogonal to the 
other parameters. Dependence on junction area A and 
capacitance Cj in Eq. (|28|l is implicit by taking the prod- 
ucts R n A = p n and R n Cj — p n tQt r /d to be constant 
in standard Josephson junctions, respectively, where p n 
is the specific (normal) barrier resistance, e r the bar- 
rier oxide dielectric constant and d the barrier thickness. 
Furthermore, we have neglected the (lnry) 1 / 3 dependence 
in Eq. (|24|l assuming ?) > 1. Hence, we can minimize 
the total circuit noise with respect to the parameters in 
Eq. i|28|) . In particular, jb appears to decrease with de- 
creasing T, a (although satisfying a ^ 1) and I c . How- 
ever, a lower I c has to be compensated by a larger L in 
order to satisfy Eq. l(T3)). for the price of lower f r cx L _1 
and a disadvantageous, although weak increase of noise 
in the second term of Eq. I|28|) . Optimum values for R s 
and R n are found from a detailed balance of the noise 
contributions. Assuming, for instance, a negligible con- 
tribution from the fourth term in J5SJ), a large R s value 
seems favorable. Due to R s <C R n , however, we see a 
conflict with a low noise requirement for the first term. 
This example implies a not too large R n j R s ratio. 



F. Dissipation and electron heating 

Without formally relating dissipation processes to elec- 
tronic noise, we consider local thermodynamics due to 
electron heating which may degrade the circuit's perfor- 
mance. The main current through the circuit is dissi- 
pated in the shunt resistor R s resulting in a permanent 
power of P s = aV s I c = R s a 2 I 2 . Successful removal of 
the excess energy is a matter of proper thermal anchoring 
(sufficiently large contact areas) to prevent overheating 
of R s . A similar power term arises during oscillations 
(in phases A and C) at the inductor Pl = V S I C . Both 
these dissipation terms essentially constitute simple heat 
loads to the cryostat. Their magnitudes are typically 
very small and the risk of local overheating of the lumped 
elements or excessive global heat overload is minor in 
standard cryogenic environments. 

However, thermal noncquilibrium in the Josephson 
junction device may have significant consequences and 
deserves closer inspection. Relevant dissipation in the 
junction can be shown to occur in phase B with Pf = 



V 2 CjV s /(2I c L) and in phase C with Pf 



VJ c /2, 
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where both are (referred to a full oscillation cycle) 
weighted with their corresponding time constants. Sat- 
isfying Eq. © yields Pf < V s I c /2, hence the heat 
in the junction is essentially P^ . In order to remove 
the hot quasiparticles from the tunnel barrier region, 
normal metal trapping layers need to be attached ad- 
jacent to the superconducting junction electrodes. At 
very low temperatures the interaction between electronic 
and phonon subsystems in a normal metal vanishes, and 
the heat flow is described by P ep = £^(T 5 — T 5 ), where 
£ w 2.4nW/(K 5 //m 3 ) is a material-dependent coupling 
constant^ v the volume of the normal metal and T e 
and T p are the temperatures of the electron and phonon 
system, respectively. The normal trapping layers will 
therefore experience an increased electron temperature 
T e w (V a I c /(2Zv) + Tp 5 ) 1 / 5 . Subsequently, the energy 
stored in the normal "cooling fins" has to be transferred 
to the substrate. However, if the temperature differ- 
ence \T e — T p \ is small and the "cooling fins" are made 
of thin films (i.e. large contact area to thickness ratio) 
the Kapitza resistance becomes negligible^ and one can 
safely assume T p to be equal to the cryostat temperature. 

IV. NUMERICAL ESTIMATES 

In order to build a relaxation oscillation circuit we are 
in principle free to choose any device or circuit compo- 
nent and derive the remaining parameters based on opti- 
mum arguments as discussed in the previous sections. As 
an example we start with a Josephson junction of given 
area A, junction (superconducting) material and an op- 
eration temperature T. The area determines R n = p n / 'A 
with a typical p n s» 1 kCl /im 2 in our standard devices 
and Cj = €Qe r A/d with e r rs 8 for A10 x and d « 2nm. 
The choice of junction material is a choice of energy gaps 
Ai, A2, determining V g , I® and Ej. It is worth noting 
that, according to Eq. I|28|) . a lower V g oc nI c R n tends to 
result in a lower noise level. A lower V g is also preferable 
to exclude perturbations like Fiske modes from the gap 
region. However, since V g delivers the oscillation ampli- 
tude, a minimum level is required for proper resolution 
of the oscillating signal V ou t ■ This conflict can eventually 
be alleviated instead by a suppression of I® by the factor 
k due to application of a magnetic field, increasing the 
oscillation frequency which may be a desired effect. Fi- 
nally, with the choice of operation point a, the required 
ratio R s /R n <C 1 and a minimum L satisfying (|I3[) . the 
full properties of the circuit are determined. 

Figure 0] shows a numerical example of a relaxation os- 
cillation circuit as a function of junction size (side length 
s), assuming an Aluminum junction (Aai = 170 /iV), an 
operation temperature T = 100 mK, the ratios R n /R s = 
10 3 , I^/Ic = 1 and h/I c — 3, and an effective L which is 
chosen 10 times larger than the minimum L mnl in (|13fl . 
The results in Fig. 21 give an idea of the order of the 
parameter ranges, including the contributions from dif- 
ferent noise sources. The flicker 1/f noise density (at 
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FIG. 4: Numerical calculations of R 3 , If> and the different 
contributions to circuit input current noise as a function of 
junction side length s. The corresponding transition current 
distributions p(i) — dP(i) /di due to thermally activated zero- 
voltage escape are shown in the bottom graph for different 
sizes s (identified by the numbers in units of pm adjacent to 
the curves.) 

/,.) was always at least 10 4 times lower than any other 
noise contribution and is therefore not shown in this ex- 
ample. The invariant parameters of this configuration 
are: f r = 5.23 MHz, uj p = 151 GHz, V s = 0.80 /xV, 
/3 C = 28.7, f3 s = (3 C (R S /R n ) 2 = 28.6 x lO" 6 , f3 L = 464.5, 
t b = 45.1 ps and R n Cj — 35.4ps. 

It is apparent in Fig. 0] that fluctuations due to ther- 
mally activated zero- voltage escape (i.e. current noise j m ) 
are the dominant noise process in the range of small 
devices. For illustration the corresponding distribution 
functions p(i) = dP(i)/di of transition currents are in- 
cluded in Fig. 01 In spite of increasing distribution width 
with decreasing junction size s, the noise density de- 
creases due to a faster decay of lb. 

As a second numerical example we calculate a realis- 
tic minimum-noise circuit configuration without leaving 
the range of classical dynamics of Josephson junctions as 
assumed for our model. We choose an A1/A10 X /A1 junc- 
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TABLE I: List of component and operation parameter values of a relaxation oscillation circuit optimized with respect to 
minimum input current noise. Initial conditions are described in the main text. 



A 


R„ 


R s Cj 


h 


h 


v s 


L 


Pc 


0L 


fr 


u> p 


T B 


<5t, 


WI 


jttt 3m 3 s 3z 


(A«n 2 ) 


(Q) 


(Q) (fF) 


(nA) 


(nA) 


(/AO 


(mh) 






(MHz) 


(GHz) 


(ps) 


(ns) 




(fA/VHS)" 


0.1 


10 4 


10 2 3.54 


25.0 


75.1 


7.51 


2.87 


26.9 


218 


104.5 


146.6 


44.9 


0.63 


0.56 


0.098 0.002*' 965 235 363 



"These units refer to all four current noise density terms. 
'Value taken at /,. . 



tion where we expect the best quality tunnel barriers and 
a sufficient oscillation output amplitude V ut ~ 340 /iV. 
The minimum junction size is restricted by the range of 
validity of our model, requiring Ej > Ec = e 2 /2Cj and 
Ej > kT to prevent single-electron charging or macro- 
scopic quantum tunneling effects, respectively* 3 ^* This is 
just satisfied with a junction of area A = 0.1 /xm 2 
(316 nm) 2 , which can be fabricated by standard e-beam 
lithography. An operation temperature of T — 100 mK is 
easily reached and maintained in modern cryostats even 
in the case of some moderately low dissipation in the 
circuit. The values for the circuit components follow 
from our definitions of a = 3, R n /R s = 100, k = 1, 
L/L m [ n = 5, and are listed in Table [I] The results show 
a total input current noise as low as about 1 pA/v^Hz, 
with the dominant contribution from thermally activated 
zero- voltage escape. Flicker 1// noise density at f r is 
at a negligible level and remains insignificant down to 
very low readout bandwidths fb w . The total noise figure 
of this configuration is well competitive with the best 
commercial SQUID amplifiers. In addition, due to the 
advantage of improving noise behavior with increasing 
oscillation frequency, it delivers a bandwidth superior to 
most SQUID systems. 

It is clear that the operation point for the current or 
voltage biased device under test is fixed to or V s in 
this exam ple. For devices requiring different bias values 
(as in Ref. 134]) the circuit components have to be adapted 
with respect to the specifications. 

As discussed in Sec. IIII Fl the total dissipated power 
P tot ps V s I c (a + 3/2) in the circuit example above (Ta- 
ble[I} amounts to about 0.87 pW. We assume the normal 
metal "cooling fins" attached to the junction electrodes 
to be thin films of 100 nm thickness and with an area of 
10x10 //m 2 . The effective electron temperature in these 
normal metal traps will then increase (from 100 mK) to 
about 108 mK, which is absolutely acceptable with re- 
spect to power equilibration in the system as well as to 
the thermal properties and operation of the Josephson 
junction. Larger tunnel junctions will require propor- 
tionally enlarged "cooling fin" areas, which is feasible up 
to a practical size compared to the bulky dimensions of 
L (and eventually R s ). Furthermore, if the normal metal 
traps also form the connecting leads, the trapping and 
cooling efficiencies may improve appreciably. 



V. POSSIBLE APPLICATIONS 

We have developed the relaxation oscillation analog-to- 
frequency converter primarily for readout of cryogenic ra- 
diation detectors^ The aim was to overcome problems or 
limitations in scaling to large number pixel readout. Be- 
sides the outstanding noise properties, a particularly nice 
feature of the relaxation oscillation circuit is its poten- 
tial for simple implementation into a frequency-domain 
multiplexing scheme by tuning the oscillation frequen- 
cies of the individual analog-to-frequency converters to 
well separated values, and then using one single line to 
read them out. It should be taken into account, how- 
ever, that a signal excursion from a detector generates a 
frequency shift, which should not overlap with a neigh- 
boring oscillator in the simplest case. A more sophisti- 
cated scheme could lock into the "dark" (no detector sig- 
nal) characteristic frequencies and, upon disappearance 
of one channel due to an analog detector signal, remove 
the other frequency bands in order to recover the signal 
of interest ~ The circuit example in Table [I] is apparently 
not a good choice for a multiplexing readout due to fairly 
broad SI m /I c . However, if we accept a moderate increase 
in noise level by choosing larger junctions, fluctuations 
are easily reduced to 8I m /I c < 0.01 (see Fig.QJ. In that 
case, and taking into account that the circuits should all 
operate in a bias range V s < V g /3, we estimate that up to 
about 10 oscillators can be implemented with sufficient 
separation. The different frequencies can be tuned by 
either different shunt resistors i? s , different inductors L 
or a careful variation of junction sizes (resulting in differ- 
ent I c ), whatever best meets experimental conditions and 
technical possibilities. One should also understand that 
combination of several circuit outputs by simple connec- 
tion (e.g. through resistors, no amplification) reduces the 
amplitudes of the individual signals by a factor equal to 
the number of interconnecting channels. Required mini- 
mum signal-to-noise therefore limits the feasible number 
of multiplexed channels. 

To test and demonstrate the working principle of a 
(single) relaxation oscillation circuit readout we have 
measured the response of a SINIS microcalorimetes2& to 
irradiation with 6 keV X-rays. The detector which re- 
placed Rb was voltage biased. Figure |S] shows the re- 
sults of an X-ray event. The circuit and device pa- 
rameters were: L = 48 nH, R s = 91m57, junction size 
,s = 15 /im and effective critical current I c — 7.28 /iA 
(k = 8). The detector's "dark" (or bias) current was J° = 
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17.5 /kA = 2.4J C , the measured analog signal peak current 
was It — 46.4 /xA, as shown in Fig. |SJi. The relaxation 
oscillation frequencies from Fig. OJ, taken at the same 
operation point and conditions, were /° = 4.47 MHz and 
fr = 12.1 MHz, respectively. Taking into account the 
conversion factor dfr/dh = 260kHz//iA, see Eq. 0, the 
analog and the frequency-modulated signal are perfectly 
compatible quantitatively as well as qualitatively (pulse 
shape). The noise level is about the same in both cases 
and is due to detector noise. The circuit noise alone is es- 
timated to contribute about 0.3 /iA rms integrated over 
full bandwidth up to f r . We should say that the mi- 
crocalorimctcr device and circuit configuration are by no 
means optimal in this example, they are rather a prelim- 
inary choice of available components. Primarily, these 
results are of illustrative nature, demonstrating the prin- 
ciple and feasibility of cryogenic detector readout. 

Other possible applications for relaxation oscillation 
based analog-to- frequency conversion can be, in a wider 
sense, considered for any type of cryogenic device oper- 
ated at relatively low bias levels, exhibiting small vari- 
ations of its electronic properties or actively deliver- 
ing small analog signals. A list may include quantum 
dots and wires, single-electron devices and quantum Hall 
structures, to name a few. Due to the large bandwidth, 
the readout method is also attractive for detection of fast 
processes like quantum noise or background charge fluc- 
tuations. The resistors Rb and R s in our scheme just 
represent a current and a voltage source, respectively, 
and can be replaced by the device of choice. 

It is important to note, however, that the oscillator 
junction characteristics (essentially represented by I c ) 
may slightly vary from cooldown to cooldown and there- 
fore cause a measurable spread in the conversion fac- 
tor Q. Therefore, the circuit is unfortunately not appro- 
priate for absolute measurements on a level as required 
e.g. by metrologists. 

As a concluding experiment we propose a setup for 
high-precision thermometry at low temperatures, replac- 
ing the classical four-point measurement on thermistors. 
The temperature-sensitive element would typically re- 
place R s to minimize dissipation. The difficulty of ap- 
plying small analog excitations and detecting low output 
levels (across long wires), competing with noise, is cir- 
cumvented by directly "digitizing" the small signal very 
close to the sensor with a low- noise converter. It is clear 
that this thermometer readout can only be operated in 
a limited temperature range (presumably one order of 
magnitude) where the junction dynamics (fluctuations, 
switching probabilities) are sufficiently insensitive to T 
variations. 



VI. CONCLUSIONS 

We have investigated the feasibility of a cryogenic low- 
noise analog-to-frequency converter with acceptable lin- 
earity over a broad range of circuit and operation param- 
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FIG. 5: (a) Response of the relaxation oscillation circuit 
modulated due to a 6keV X-ray event in a SINIS detector 
at 80 mK. The sinusoidal shape and reduced amplitude of 
the oscillations are due to band-pass filtering at the post- 
amplification stage, (b) Rectified oscillator signal where each 
line represents the position of an oscillation cycle maximum 
(note the larger time scale, applying to all three lower plots), 
(c) Inverse oscillation periods as a function of time, equiva- 
lent to a time-dependent frequency f r . (d) Analog signal from 
the same detector measured with conventional electronics and 
taken at corresponding experimental conditions. 



eters. The dynamical behavior can be well described by 
simple circuit theory and classical models of the single 
Josephson junction involved. Their agreement with ex- 
perimental data is perfect. We have presented a thorough 
analysis of noise sources, where fluctuation processes in 
the Josephson junction appear to usually dominate the 
circuit's noise figure for typical configurations and exper- 
imental conditions. The inherent broadband operation 
paired with very good noise performance offers a ver- 
satile system for a wide range of applications. As one 
possible example we have demonstrated the readout of 
a cryogenic microcalorimeter exposed to X-rays. Imple- 
mentation into a multiplexing scheme was discussed and 
needs to be experimentally tested for a large channel- 
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number readout. 
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suppressed. Setting h = 0.5 in Eq. (jA.lfl yields the ex- 
pectation value 



< 7 ™> ~ i 1 fo i 
— - — = mi « 1 27 In , , 

Ic 8 V ln(9?/); 



2/3 



In the limit (77 — > 00) this expression approaches 
%l«l-i(2 7 ln, ? ) 2/3 . 



(A.2) 



APPENDIX 

For the analysis of noise due to thermally activated 
zero- voltage escape an evaluation of the central moments 
of the probability density function p(i) — dP(i)/di is re- 
quired, yielding expectation value mi = fip(i)di and 
variance 7712 = /(* — mi) 2 p(i) di, where P(i) is given by 
Eq. I|22|l and a m = ^/ma subsequently denotes the stan- 
dard deviation. Analytical integration can be circum- 
vented, however, by approximating p(i) with a Gaussian 
distribution and solving for the appropriate values satis- 
fying P(mi) = 0.5 and P{ mi ±a m ) = 0.5{l±erf(l/ v / 2)}, 
respectively. A general solution of P(m) = h is 



m 



2jW 



ln(l - hf 



2/3 



(A.l) 



where r\ — 2(uj p /2irf r ) 6 ('y/2) 5 and W is the Lambert 
function satisfying W(x) ■ exp(VF(x)) = x. An asymp- 
totic expansion of W for (x — > 00) can be written as 
W(x) — In a; — ln(lnx), where higher order terms are 



Correspondingly, for (6I m ) jl c = a m we find an expan- 



4.73 1.58 In (In 77) - 10.15 



i(2 7 ) 2/3 (InrjjVS (ln?7) 4 /3 
which, for (77 — > 00), approaches 

(Sim) r l r2V> a/3 4 -73 
Ic 8 1 11 (ln 7? ) 1 /3 



0.94- 
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(A.3) 



"(lnry) 1 ^ ' 

If m\ significantly deviates from unity the approximation 

describes the relative distribution width of transition cur- 
rent fluctuations. 
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